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Abstract. Recent Hcrschel/SPIRE maps of the Small and Large Magellanic Clouds 
(SMC, LMC) exhibit in each thousands of clouds. Observed at 250 microns, they must 
be cold, T ~ 15 K, hence the name "Herschel cold clouds" (HCCs). From the observed 
rotational velocity profile of the Galaxy and the assumption of spherical symmetry, its 
mass density is modeled in a form close to that of an isothermal sphere. If the HCCs 
constitute a certain fraction of it, their angular size distribution has a specified shape. 
A fit to the data deduced from the SMC/LMC maps supports this and yields 1.7 pc for 
their average radius. There are so many HCCs that they will make up all the missing 
Halo mass density if there is spherical symmetry and their average mass is of order 
1O,OOOM . This compares with the Jeans mass of circa 40, 000M Q and puts forward 
that the HCCs are in fact Jeans clusters, constituting all the Galactic dark matter and 
much of its missing baryons, a conclusion deduced before from a different field of the 
sky (Nieuwenhuizen, Schild and Gibson 2011). A preliminary analysis of the intensities 
yields that the Jeans clusters themselves may consist of some billion MACHOs of a 
few dozen Earth masses. With a size of dozens of solar radii, they would mostly not 
lens but cause occultation of stars in the LMC, SMC and towards the Galactic center, 
and may thus have been overlooked in microlensing. 

PACS numbers: 

95.35.+d Dark matter (stellar, interstellar, galactic, and cosmological) 

98.35.Gi Galactic halo 

98.38.Jw Infrared emission 

98.56.-p Local group; Magellanic Clouds 
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1. Introduction 

The stone which the builders rejected, 
the same is become the head of the corner 

Matthew, XXI: 42 

The dark matter riddle was formulated for the Galaxy in 1922 by Jacobus 
Kapteyn [1]: there must be more mass than in stars, the invisible or dark matter. 
Currently it is mostly approached from particle searches, aiming to find the weakly 
interacting massive particles (WIMPs), the purported cold dark matter (CDM) particle 
[2]. A complementary approach is to investigate statistical properties of dark matter. 
Since the Galaxy exists billions of years, its dark matter, whatever its nature, probably 
had time enough to reach some type of equilibrium. This expectation is supported by 
the observed flattening of galactic rotation curves. After all, a really constant circular 
rotation velocity v c would imply, due to the Kepler relation v c = [GM(r)/r ] 1//2 , that 
the mass M{r) inside a sphere of radius r grows linearly in r, and hence that the mass 
density decays as r~ 2 . But that is the well known case of the singular isothermal sphere, 
a solution to the Poisson-Boltzmann equation, where the mass density is a Boltzmann 
distribution, the equilibrium distribution of standard statistical physics. The Galaxy is 
close to that. 

The question "Which type of dark matter causes the isothermal profile?" has no 
established solution at present. But galaxy-scale difficulties for WIMP dark matter are 
abundant, see, e. g., [3]. It is proper to approach the problem in the opposite direction: 
identify a candidate, and investigate whether it satisfies the relevant statistical criteria. 

The answer to this may be related to the missing baryon problem: of all baryons 
expected from the standard model, 10% are known to be in stars and 20-30% in gas of 
galaxies and galaxy clusters, while the rest, 60-70%, is still unaccounted for [4]. A new 
component, hot gas in the Galactic halo that weighs over 10 billion solar masses, has 
been reported recently [5], but it will likely not be enough to solve the full problem. 

Our argument starts with the remarkably low 15 K temperature of clouds all over 
the sky, known from IRAS observations, and commonly attributed to dust, but see [6]. 
For lack of suitable instruments, 15 K objects would have mostly remained in the dark 
until recently, but by now Planck detects 915 cold molecular cloud core candidates, out of 
over 15,000 unique sources that litter the sky at the 857 GHz (350/x), the largest Planck 
frequency, see Fig. 6 of [7]. Several of these clouds have a gas mass of 10,000 solar masses 
[8, 9]. Towards the Galactic center, at coordinates \b\ < 1° and i = 300° - 330°, 1205 
dark clouds are detected by Herschel/Spitzer, which can be seen not only in absorption, 
but also in emission at 250/i and 500/i [10]. 

One particular 15° x 12° region on the sky observed already with BOOMERANG 
in the 1990's and focused on recently by Veneziani et al (2010) [11], exposes dozens 
of 7-20 K clouds. Their temperature, typical opening angle, two-dimensional number 
density and emission intensity are explained from gravitational hydrodynamics (GHD) 
employing an isothermal sphere of Jeans clusters (JCs) in which, next to stars and 
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Figure 1. This luminosity- inverted 8.32° x 4.16° Herschel-SPIRE map of the SMC 
at 250 /Ltm exposes the star-rich bar on the right and the wing towards the left, and 
thousands of Herschel cold clouds (HCCs) in the Galactic foreground, of different size 
and intensity. This figure was prepared from data made public by NASA/ESA. 



gas, all dynamical mass is contained [12]. It was concluded that these clouds are 
JCs and make up the galactic dark matter and most of the missing baryons. This 
challenges the standard ACDM or concordance model, and clearly awaits testing by 
further observations. 

Cold 15 K objects should be clearly exposed by Herschel SPIRE's 250/i detector. 
In early 2012 NASA/ESA published a multi- wavelength map of the Small Magellanic 
Clouds (SMC) [13]. (A similar picture for the Large Magellanic Cloud [14] is less suited 
for our purpose; it will be discussed in section 5.3.) In the SMC picture one observes 
the bar on the right and the wing from the right to the middle, which both contain 
many stars. Our interest is in the 250/i maps, represented as red dots in these color 
figures, that we reproduce for the SMC as gray dots in Figure I. Surprisingly, one 
observes many hundreds, if not thousands, of HCCs all over the picture (except in 
regions without data), not exposed before in this number. Since most of them seem 
absent in the "green" map (160 and 100 /iin), they are likely not background galaxies, 
but must be cold (~ 15 K), and we adopt the Juvela (2012) name "Herschel cold clouds" 
(HCCs). They are the focus of our present study. Our aim is to investigate whether the 
HCCs can correspond to JCs that would comprise the full dynamic dark matter of the 
Galaxy [12]. 

In section 2 we discuss the Jeans theory and its extension by gravitational 
hydrodynamics, and we recall a number of observations already explained by this model. 
In section 3 we formulate a model for the expected angular size distribution of the HCCs. 
In section 4 we investigate their size distribution in the SMC map and in section 5 we 
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discuss the fit to our model. In section 6 we discuss the results within the view point of 
JCs and micro brown dwarfs (/iBDs), and we conclude with a discussion. 

2. The Jeans instability and GHD 

The Jeans 1902 mechanism asserts that when in the early Universe the plasma changes 
into a H-He gas, this gas, considered to be infinitely distributed, exhibits instabilities 
at the Jeans scale, where pressure can no longer overcome the effects of gravitation and 
diminished heat transfer [15, 16]. This leads to the formation of Jeans gas clumps that 
weigh about 40,000 M . Globular star clusters are supposed to have formed from a 
handful of coagulated clumps, but other ones are generally considered to have turned 
into stars, brown dwarfs, planets, gas clouds, etc., so that they should not exist anymore, 
at least not in large numbers. We shall pose, however, the question: does the majority 
of the Jeans gas clumps still exist as HCCs? 

In GHD it is asserted that, due to turbulence and viscosity, the Jeans clumps 
themselves fragment in earth mass MACHOs, called primordial gas particles [16] or, 
what we shall employ here, micro brown dwarfs (/iBDs) [12]. This turned the Jeans gas 
clumps into JCs of billions of /iBDs. In order to reach a temperature of 2.725 K, the 
/iBDs would have to cool from gas to solid, by radiating away the latent heat in order 
to cross the first-order phase transition line that emerges the H triple point at 13.8 K. If 
this takes more than a Hubble time, their temperatures must become stuck somewhere 
above 14 K. 

The wide-binaries problem, the last-parsec problem and the Helium-3 problem also 
find a simple explanation in GHD [17], as do long-duration radio transients [12, 18, 19]. 
GHD offers plenty of dark baryons for making new stars even with modest recycling 
of stellar matter [20, 21, 22], because it would consider the outflowing and infalling 
matter as dynamical behavior of the galaxy's large baryonic halo. The /iBDs would be 
leaky, so there is gas throughout the JC. This offers an explanation for the Lyman-alpha 
forest: thousands of extinction lines in quasar spectra due to hydrogen clouds along the 
sightline. Why these clouds did not dissolve quickly is still a miracle [23], but JCs of 
leaky /iBDs provide a simple explanation: the gas is bound by mass concentrations, the 
/iBDs. Therefore, the total mass can be much larger than the gas mass, and this is the 
picture that we shall investigate in the present paper. 

The idea can also be tested in our immediate neighborhood. The Local Leo Cold 
Cloud is a large (~ 5 deg), cold (15-20 K), nearby (10-20 parsec) H cloud [24]. Like the 
Lyman-alpha clouds, it can be long lived if it contains many unresolved /iBDs, which 
bind the gas gravitationally, thus explaining why the gas can have the observed large 
pressure discussed in [25]. The observed velocity differences inside the Local Leo Cold 
Cloud are explained by the local winds due to mostly random motion of /iBDs, that 
sweep the gas with them. 

In GHD, planet and star formation occurs by coalescence of /iBDs [12, 26, 27]. 
They are cold and relatively small now (solar size), but when heated they expand and 
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can coagulate. First a gas planet is formed and, after enough steps, a star is created. 
This explains why galaxy mergers such as Antennae [28] and Tadpole [29] have tails 
of thousands of young (< 10 Myr) globular star clusters, much younger than the event 
duration of galaxy merging (150-200 Myr): during the merging, the local JCs of one 
galaxy are heated by the center of the other, so that they later develop stars [12]. Many 
or possibly all stars may have been formed in such clusters [28]. 

2.1. The Jeans mass 

We shall need a mechanism that explains the formation of so many, so regular, so 
identical cold clouds. A well known candidate is the 1902 Jeans instability of the 
neutral gas, just as it forms at the recombination. It is often taught as a separate 
issue in standard courses on astrophysics, but for us it will be a cornerstone. 

The Jeans mechanism asserts that at last scattering (L; decoupling, recombination) 
the plasma turns into a neutral gas, which is unstable on a scale where pressure and 
heat transport are overcome by gravitation. The free fall time is r g = (Gps)^ 1 ^ 2 = 5.30 
Myr, where the baryonic density at the decoupling at zl = 1088 is p B = ^bPc(1 + -2l) 3 — 
5.35 10~ 22 gr cm -3 for Q B = 0.045 and p c = 9.21 10 -30 gr cm -3 at h 70 = 1. The speed 
of sound of a monoatomic gas is v s = (5p/3pY^ 2 . For the gas of H and 25% of weight 
in He the pressure p = 0.8l3pk.BT/m N yields v s = 5.68 km s -1 . The gas fragments at 
the Jeans scale Lj eans = v s r g = 31 pc into Jeans clumps with mass 

M Jcans = -paeans = g ^Pc(l + ^)^Lns = 6g3/2 ' 1/2 = 40 ' 000 M «" (*) 

A practical way to check this is to consider globular star clusters, which, we presume, 
are the result of up to a dozen or so of JC mergers, as is exhibited by their sets of stars 
of different ages [30]. To yield a GC mass between ten thousand and a few hundred 
thousand solar masses, the Jeans mass should indeed be a few ten thousand solar masses. 

2.2. Nearly isothermal mass density profile 

The circular rotation speed at the position of the Sun, taken as v* = 220 km s _1 , drops 
within 10% accuracy mildly to 175 km s _1 at 60 kpc [31]. This range can be modeled as 
v c (r) = v*(d SA /r) £ km s" 1 with e = 0.114 and d SA = 8.0±0.5 kpc, the distance to SagA*. 
Now it is found that v* = 240 ± 14 km s" 1 [32]. In the presence of spherical symmetry, 
Kepler's law GM(r)r _1 = f 2 (r) then yields for the mass density p = M' {r) / Anr 2 : 

^ =il ^{¥) 2 '- < = 240 ± 14km S ->. (2) 
For e = this is the singular isothermal sphere model with velocity dispersion v*/\/2. 
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3. Model for the angular size distribution of Galactic foregrounds 

We assume that the HCCs make up a fixed fraction Ch cc < 1 of the total mass density. 
With average size R and average mass M, this leads us to their average number density 

n[r) - u p[r) _ Rr2 y r j , A- . {6} 

For the distribution of radii R at position r we consider the factorized form n(r, R) = 
n(r)f(R/R)/R with / °° dx f(x) = J °° dxxf(x) = 1, which reproduces the local average 
/ dRn(r, R) = n(r), while evidently (R) = R. We specify to the T-distributions f(x) = 
(nx) n exp(— nx)/xT(n), which for n — > oo lead to the fixed- R case f(x) = S(x — 1). 

3.1. Seen from Sag A* 

From SagA*, an HCC of physical radius R is seen with angular radius 9 = arctan R/r m 
R/r and angular diameter 29. Let N(9) denote the number of HCCs that have an 
angular radius smaller than 9. Its differential increase per unit angular radius per unit 
solid angle is 

where dfl is the differential solid angle, the integral of which cancels the factor AQ. 
Using 5(9 — Rr" 1 ) = R9~ 2 5(r — R9~ l ) and defining x = R/R, the result takes the form 

diV(#) AC (6) [™ . , fd SA 9\ 2£ r(n + l-2e) fnd SA 9 x2 



d9dft 9 2 



C(9)=j o dx/(x)x(_j = i — \ ,(5) 



which has dimension rad -3 . Notice that A and R do not occur independently here. 
For e — > one considers in fact isothermally distributed HCCs, where this relation has 
C — 1 and expresses the isothermal relation p(r) ~ r -2 in a dimensionless form. 

3. 2. Seen from the earth 

The fact that the Sun is not located in the center of the Galaxy introduces a slight 
complication. But actually it turns out as a benefit that allows us to estimate both the 
amplitude A and the average radius R from the statistics of the data. 

A sphere of radius R = xR at distance d from the Earth at galactic coordinates 
(6, £) has position r = (d cos b cos £ + dsA , d cos b sin £, d sin b) with respect to the center of 
the Galaxy. The modulus is r(d, b, £) = (d 2 + 2d§ A d cos b cos £ + g^a) 1 ^ 2 and the volume 
element d 3 r = dddVld 2 , where dVt = dbd£ cos b. Because 9 = R/d we have the definition 



dN{9) 

d9dtt ' Ail 



1 —JJdddQdxd 2 n{r)f{x)5(9- X ^\. (6) 
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When AQ is a small field around (b, £), the integrand hardly varies with Q, so AQ drops 
out again. The d-integral yields dN(9)/d9d£l = A9~ 2 C(9), with A from (3) and now 



4. Data analysis of Herschel maps of the Magellanic Clouds 

4-1. Detection 

The Herschel "Spectral and Photometric Imaging Receiver" (SPIRE) instrument has 
detectors for central wavelengths of 250fi, 3G0fi and 520fi. The full-width at half- 
maximum of the 250/i detector is 18.1", see section 6.4 of Ref. [33]. Herschel-SPIRE 
has delivered a multi-frequency map of the SMC [13]. The 250/i map has a resolution of 
6000 x 3000 pixels and its angular size is determined by identifying a number of spots 
with stars known from catalogs. We thus obtain 8.315° x 4.157°, which corresponds to 
pixels of linear size of 9 P i x = 4.99". Hence the resolution limit is 3.6 pixels across. 

We shall present data for structures from 8 to 124 pixels in size. Each pixel has a 
solid angle 9 2 ix . Approximating an object with a number of pixels iVp ix as a sphere, we 
define the angular radius 9 by N pix 9 2 ix = 1\6 2 . In the analysis the angular widths are in 
the range 15" < 29 < 63", the smallest value being just below the resolution limit, but 
not causing a clear deviation in Fig. 2. 

4-2. On Image J 

To count the HCCs we use the the freely available program "Image J". We first split the 
colors of the multi-frequency map to extract the 250/i map represented in Fig. 1. Next 
we mention some relevant points/decisions made with respect to the counting method. 

1) First of all, we set the minimum threshold value to 11 and the maximum to 255 
(on a 256 bit gray-scale) to select the pixels with a grayscale between 11 and 255. If the 
minimum threshold value is set much higher, visible absence of data will occur; setting 
the threshold much lower than 11 causes too much background noise. Here we define 
background noise as very faint objects measuring less than 2x2 pixels, which is below 
the diffraction limit of ~ 4 x 4 pixels for Herschel's 3.5 meter primary mirror at this 
wavelength. 

2) We use a watershed function to split overlapping HCCs. This also helps the edge 
tracker (see point 5) to distinguish between the background and HCCs. 

3) We count the number AiV of HCCs for pixel sizes between 6 and 125. 

4) Large distortions (defined as HCCs bigger than ~ 130 pixels), like luminous 
matter clumps, are unusable. They occupy 45-55% of the map, while some 12% of the 
map is without any data. In such regions we normalize by the analyzed area, see 6). 

5) For the counting we use an edge tracker to help distinguishing near by /overlapping 
HCCs. However, it seems that the edge tracker increases the apparent area size of HCCs 




(7) 
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as a function of their true area size. We correct for this issue in our data by the empirical 
relation N^ c = (0.11 ± 0.04) (A^ 61 *) 1 - 41 * - 01 . 

6) We bin the data in sets of AN pix = 4 pixels, viz. iV pix = (6, 7, 8, 9), (10, 11, 12, 13), 
etc., with central bin values iV p i x = 7.5, 11.5 etc., to be used in the definition of 9. We 
determine 6 2 &N/&6dVl = 2ir- 1 / 2 N^6 pix ■ AN/AN pix Att, where AN is the number of 
HCCs in the bin of the investigated region of angular area Af2 = cos bAbAi. 

7) We determine the error bars by measuring 7 different areas in our data set and 
calculating the standard deviation for each bin of 4 pixels. 

4-3. Two ways of analysis for the SMC 

We divide the SMC map into 7 "good" regions that are poor in stars and define region 8 
as the rest, centering on the star-rich regions, on the border regions without data, and on 
other regions with poor data. The average angular position is (b,l) = (—44.2°, 302.8°). 

In our first method, we analyze the full SMC map. ImageJ treats the luminous 
regions as big HCCs of 10,000 pixels and more, so they do not add to our counts. We 
correct the useful area by subtracting these big luminous regions, as well as regions of 
HCCs with more than 130 pixels. We do the same for the side regions without data. The 
error in the variable y = 8 2 dN(9)/d8dfl is defined as Ay = [X)™ =1 Wi(yi — y) 2 Yn=i w, 2 ] 1 ' 2 ) 
where n = 8 and Wi = AQi/ Z)" =1 AQj is the areal weight of the region % and y = WiUi- 

In our second method for the SMC, we select the seven "good" areas away from 
the star-rich regions, we analyze them and we join them by area size. The error in the 
variable y then sums over these n = 7 areas. 

For the LMC only a strip on the left and a smaller strip on the right are useful for 
our analysis. The errors are now larger, since we only have two regions (n = 2) with 
relatively few HCCs. The average position of the LMC is (b,E) = (-32.89°, 280.47°). 

4-4- Provisos of our method 

We take for granted that the "red dots" of the color version of the SMC map (the gray 
ones in Fig. 1) represent physical objects and we model them as spheres. One may argue 
that the observations present more edgy or hairy objects, which could be an effect of a 
mutual interactions, near passages or central collisions, or complete absence of nearly 
spherical objects. But most of these issues concern details below the resolution limit, 
which cannot be analyzed properly. In the case of large objects, we assume that they 
are agglomerations of spheres, and separate them accordingly. But there can be mass 
streams between them, in particular due to mutual interactions. However, given the size 
of our statistical errors, these issues do not seem too relevant at our level of description. 

The LMC map [14] is largely dominated by the star-rich region. For our analysis 
small side regions remain available, the results of which agree with those for the SMC. 

Although we consider our model as reasonable and the parameters come out as 
reasonable, we stress that the data could be fit to many other models. 

We have analyzed only the MC regions, which cover a small piece of the sky. 
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o-^bd(km/s) 


SMC-total 
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23.6 


0.84 


102±4 


1.68±0.18 
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4.16±0.07 


SMC-good 


oo 


67.1 


2.40 


106±3 


1.82±0.18 


6.2±0.6 


14.1±1.8 


4.08±0.06 


LMC-sides 


oo 


12.5 


0.45 


113±14 


1.17±0.19 


5.0±0.8 


8.6±2.4 


3.97±0.25 


S+L MC 


oo 


38.1 


0.66 


i r\ A i o 

104±3 


1.65±0.15 


5.6±0.5 


13.1±1.6 


4.14±0.07 


SMC-total 


2 


22.4 


0.80 


110±5 


1.30±0.18 


4.4±0.6 


9.8±1.7 


4.02±0.08 


SMC-good 


2 


65.5 


2.34 


114±4 


1.43±0.18 


4.8±0.6 


10.4±1.7 


3.95±0.07 


LMC-sides 


2 


13.4 


0.48 


127±20 


0.79±0.18 


3.4±0.8 


5.2±1.9 


3.74±0.29 


S+L MC 


2 


37.5 


0.65 


112±4 


1.25±0.14 


4.2±0.5 


9.2±1.4 


3.98±0.08 



Table 1. Fits of the model for Jeans clusters for the four data sets: the total SMC, 
the good regions of the SMC, the sides of the LMC, and the total SMC combined with 
the sides of the LMC. A is an amplitude; R the average radius of the HCC; #mc refers 
to #smCi the angular radius of an average HCC located at the SMC, except for the 
LMC-sides cases, where it correspondingly refers to the LMC. Finally, M is the average 
HCC mass needed to explain all dynamical mass in the Galaxy and c^bd is the velocity 
dispersion of the yuBDs. The upper data refer to HCCs with a fixed radius (n = oo) 
and the lower to radii taken from the n = 2 T-distribution f(x) = 4xcxp(— 2x). 

In this study we did not have access to the infrared intensities, which give additional 
information about the ratio of their radiating surface to their distance squared. 

5. Fits to the nearly isothermal model 

5.1. Towards the SMC, the case of equal radii 

A x 2 fit of the theory of previous section to the observed 30 binned values of 6 2 dN/ d6dQ 
has been performed, see Fig. 2. For the full SMC picture the fit yields 

y 2 

A = 102 ±4, R= 1.68 ±0.18 pc, x = 23.6, ^- = 0.84 (8) 

where ^ = 30 — 2 is the number of degrees of freedom. The error bars derive from 
the inverse correlation matrix of the fit parameters p± = A and p2 = -Rpc -1 for the 
prediction = AC(9i;R), viz. = 2Y,i(dyi/dp a )(dyi/dpb)o'^ 2 . It follows that 

(5A 2 ) = C n = 22.0, (5A5R)pc~ 1 = C 12 = -0.71 and (5R 2 )pc' 2 = C 22 = 0.032. 

The angular radius of a HCC of average size, located in the SMC, is ^smc = R/dsMC, 
with g?smc = 61 ±3 kpc the SMC distance. Hence ^smc = (5.7±0.6)", which is indicated 
in Fig. 2. (By linear regression, we express here and below all errors in terms of C a b-) 

To test whether the HCCs can belong mainly to the SMC, we follow our second 
method and discard the data from the star-rich SMC region, about 51% of all data. But 
this appears to yield only minor shifts in A and R, see the second row in Table 1. 
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Figure 2. N(6) is the number of Herschel cold clouds (HCCs) with angular radius 
smaller than 9. The data for 9 2 dN/d9dfl (in rad -1 ) as function of 9 (in arcsec) are 
deduced from Fig. 1. The full curve is the best fit for HCCs of equal radius and 
the dashed curve for HCC radii taken from a T-distribution with n = 2. The initial 
growth stems from the 9 2s factor and the decay from the power of dr -1 in Eq. (7). 
The vertical bar denotes the angular radius 5.7" of an HCC of average radius located 
in the SMC. 



Both the small values of #smc an d the minor shifts with respect to leaving out 
star-rich regions reject the assumption that the HCCs of Fig. 1 belong to the SMC. 
Hence we conclude that they nearly all belong to the Galactic Halo. 

The investigated distances stem from R/8 and lie between 10 and 40 kpc. All fit 
parameters are tabulated in Table 1. As a canonical average radius, we adopt R = 2 pc. 

5.2. Towards the SMC, spread in radii 

The radii of the HCCs can be random, mixing angular sizes of large, far HCCs 
with nearby, small ones. As an example, we consider the n = 2 T-distribution f(x) = 
4xexp(— 2x). Rows 5, 6 and 8 in Table 1 show differences with respect to the equal- R 
cases (rows 1,2 and 4, respectively). These modest deviations make it unlikely that the 
size distribution can play an important role in the physical interpretation of the results. 

5.3. Towards the LMC 

A similar map exists for the LMC [14] and in the data acquisition, described above, we 
only use two side regions. Now the error bars are much larger, which actually leads to 
a small \ 2 (Table 1, rows 3 and 7). We consider the deviations from the fit parameters 
for the SMC (Table 1, rows 1, 2 and 5, 6) to be statistically unimportant. 
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5.4. SMC combined with LMC 

Combining the 30 data for the total SMC with the 30 for the LMC, hardly changes the 
parameters of the total SMC alone and decreases only slightly the error bars (Compare 
rows 5 and 8 in table 1 with rows 1 and 4, respectively). This shows again that the 
LMC data, with their larger error bars, have low statistical weight. 

6. Interpretation in terms of Jeans clusters and micro brown dwarfs 

Recalling Eq. (3), we get for the average mass of the HCCs from the SMC 



The first (second) estimate holds the assumption of fixed radius for the full (restricted) 
data set; for the case of random R this becomes (9.8±1.7)1O 3 M or (10.4±1.7) x 1O 3 M 
for these sets, respectively. For the LMC the results deviate, see Table 1, but, as 
discussed, they have lower statistical weight. 

Without further input, Ch cc m ay take a value of order 10 -3 — 10~ 4 , the case of solar 
mass gas clouds, as often concluded from gas masses alone. However, if the clouds are 
composed of leaky yuBDs, the total mass can be much larger. 

The Jeans mass has been estimated as 4O,OOOM in Eq. (1). Therefore, 
it is imperative to discuss the case Ch cc = 1, for which the HCCs represent the 
complete dynamical dark matter of the Galaxy. One then identifies the HCCs of 
heretofore unknown origin with the JCs predicted by the Jeans instability and the 
GHD fragmentation into /iBDs. In view of table 1, we adopt the canonical average mass 
Mj C = 1O,OOOM . This lies within the range of HCC masses deduced from gas alone 
[8, 9, 35]. The typical distance to the nearest JC is d\ = [p(dsA)/Mj c ] _1 / 3 = 92 pc, its 
angular radius being 1.1°. Cold structures of size ~ 0.5° have been discussed in [11], 
while the HCCs analyzed in [8, 9] have distances between 225 pc and 1.67 kpc. 

If the JCs themselves are isothermal distributions of /zBDs with velocity dispersion 
o" M bd, their mass and radius are related as GM- ]C = 2cr^ hd Rj C . Assuming that <r M bd is 
basically a constant, this relation also holds for the average. We can thus express A from 
(3) with c hcc = 1 as A = (1 -2e)vf /8nal hd , which yields, for fixed-/?, a Mbd = 4.16±0.07 
km s" 1 (4.08 ± 0.06 km s" 1 ) for the first (second, restricted) SMC data set, and when 
modeled with random radii <x M bd — 4.02±0.08 km s~ x (3.95±0.07 km s -1 ), respectively. 
The LMC yields 3.97 ± 0.25 km s" 1 for fixed R and 3.74 ± 0.29 km s" 1 for random R. 
As canonical value we adopt the value a^bd — 4 km s -1 , known from stars in globular 
clusters. The value Ch cc = 1 is also desired to avoid a too small cr^bd ~ 4c^ km s^ 1 . 

6.1. On intensities 



M Rvf(l - 2e) 
c hcc A-kGA 



(13.5 ± 1.9)1O 3 M or (14.1 ± 1.8)10 3 M ( 



(9) 



At the time of this research, the intensities of the SMC map were not made public, so 
they were not considered. By now the rough data are available on the Herschel site and 
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clean data are to be produced from them (M. Meixner, private discussion). In a future 
study we intend to incorporate that information for a further test of our model. 

Here we look at a few characteristics. The Planck function gives the amount of 
electromagnetic energy radiated by a black body in thermal equilibrium per unit of 
time, frequency, area and solid angle, B U {T) = 2hu 3 c^ 2 [exp(hh / /k B T) — l] -1 . If the 
/xBD surface is optically thick, one expects a nearly exact Planck spectrum. The /zBDs 
in a given JC may have a distribution r(T) of temperatures, which leads to B u = 
J (YTt(T)B v (T). In the modeling as cosmic dust, one often takes B v « (v /vq) 13 B u (Tq) 
with spectral index (3^2 [8]. Here we continue with assuming a pure Planck spectrum. 

The radiation into the outward hemisphere picks up an angular factor tt. The energy 
emitted per second per Hz by A^bd thermal spheres of radius R^hd is 4ir 2 N^B^R 2 ^ 
Measured per unit area at distance d, this amounts to a specific intensity 



T D ^W-R^bd D ^jc-R^bd D ^bd-Rjc-R/xbd D ^bd-^jubd 2 /ln x 

'' =fC "^ =lB "v = * B " ="*-4^rs- (10) 

In the last identities we model, next to the JCs, also the /iBDs as isothermal spheres [12]. 
They are supposed to consist of gas, H with the primordial He fraction = 25% in 
weight. Neglecting metals, this leaves pu = rinm^ = (1 — (f>)p. Withp = (nn + nu c )kBT, 
we denote <r| as = p/p — (1 — |0)A;bT '/tojv, which yields cr gas = 317ms _1 (T/15K) 1 / 2 . 

In a statistical model we may define the intrinsic JC intensity i„ = I v d 2 R~ 2 = I u 0~ 2 , 
and assume local fluctuations with density f v (x,i„) normalized as / dif u (x,i) = f(x). 
This yields dA^/d^dfid/^ = A$- 4 C V (9, l v Q- 2 ) with C v (9, i) from (7) with f(x) ->■ /^(x, i). 

In the raw 250/i Herschel data each pixel has size #p ix = 6" and a saturation 
value between and a few, in units of Jy beam -1 . We denote by Sj C the average 
saturation value of the pixels that expose the considered JC. The effective beam area 
is fibeam = 9.28 1CT 9 sr [34], so the intensity recorded by a pixel with saturation Sj C is 
91.2 Sj C mJy. A JC of angular radius 9 involves a number of pixels = n9 2 9 p 2 . and an 
intensity I u = Sj C ii9 2 Q^ aia Jy. Equating this with (10) yields an estimate for the /iBD 
radius, 



_ s jc sr Jy a 2 gas s jc T ff„(15K) / 4kms 1 V R ]C 

Mbd Abeam B„ <rJ M JC °0.05 15K B U (T) { a, hd ) 2pc' [LL) 

i.e., dozens of solar radii. The /iBD mass then follows as a good dozen of Earth masses, 



With typically Sj C ~ 0.05 this compares with the estimate of 13M e [36], and is somewhat 
larger than the ~ 3M e deduced from quasar microlensing observation [37]. 

We can also consider the microlensing of an SMC star by an /iBD inside one of the 
JCs in front of the SMC, at distance d = xd SMC with x < 1. The Einstein radius reads: 
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7? - 2 frh rVUf A -n-Uff °> d 4a:(l-x) B V {T) 2pc 

" c V^ 1 " *) Gil WsMC - 0.34i? Mbd — 2QSjc ^(isK)/^ 



(13) 



(For a star towards the Milky Way the factor 0.34 would be at most 0.12.) Unless 
T is well above 15 K, one typically has Re < R^bd- This is the case of occultation, 
where the big fiBD absorbs the direct light when it passes in front of an SMC star, 
while its outer atmosphere refracts light [38]. With the /zBD speed dominated by that 
of the respective JC in which it is embedded, the lensing-occultation event would last 
t occ ~ 2R^ hd /v* ~ 3.5 days. 



7. Conclusion 



Early in 2012, an intriguing multifrequency map of the Small Magellanic Cloud was 
made public. In its color presentation, it contains thousands of red dots, Herschel cold 
clouds (HCCs), related to the 250 micron observations that we present in Fig. I. 

As a second year project at the University of Amsterdam, it was proposed to analyze 
the distribution of these HCCs using the program "ImageJ". With only the angular 
shapes available, no statement can be made about individual objects, so the approach 
has to embody a statistical analysis (discussed in section 4) and a comparison with 
theory. Preliminary results were reported at the end of the term [39]. 

In the theoretical modeling we assume that the Galactic halo is essentially 
spherically symmetric, taking into account that the dark matter, whatever its nature, 
is probably well mixed. A statistical analysis is promising, because it exhibits a scaling 
behavior which, if valid in some region, may apply to the full Halo. 

Since we count many (6600) HCCs in the 50% usable part of the ~ 8° x 4° SMC field, 
we assume that they represent a certain fraction of the total mass, and are well mixed 
too, so that their distribution has the same shape as the dynamical mass distribution. 
The Galaxy's slowly decaying circular rotation speed is modeled as v c (r) = v *(g?sa/V) £ 
with e = 0.11, which corresponds to the nearly isothermal (~ r - 2 ~ 2e ) mass density for 
the Galaxy, and, in the model, to an r~ 2 ~ 2e number density for the HCCs. We allow 
that the density has also local fluctuations in the HCC size. The number of HCCs 
per unit angular radius per unit solid angle, dN(9)/d9dil, behaves roughly as A9~ 2 for 
5" < 9 < 30". The amplitude A ~ 100 is reasonable, taking into account that typically 
9 ~ 20". This corresponds to an amount of HCCs in the field of Fig. 1 of the order 
8° x 4° x dN/dQ ~ 0.010 dN/d9dtt ~ 0.01 A/ 6 ~ 10, 000, indeed about twice the number 
6600 of HCCs with N pix > 6 that we count up to 40 kpc in the usable areas. 

Another intuitive interpretation of our numbers is as follows. If we take our SMC 
count of 6600 HCCs, multiply this by 5 for covering the Halo up to 200 kpc and then 
by 4tt/ (50% x 8.3° x 4.1°) = 2400 for covering the whole sky, we arrive at an expected 
total number of 7.9 10 7 HCCs for the Galaxy; if they weigh 10, 000M o on average, 
this constitutes 0.79 10 12 M Q in total. This rough estimate stems with the expected 
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total Galactic mass, and can be compared with the 1.11O 12 M from Eq. (2), that is 
incorporated in our detailed modeling of the spatial HCC distribution. 

Seen from the Earth, nearby objects look larger than from SagA*, which induces a 
depletion of the HCC histogram towards larger angular radii, up to a factor 2 at 30". 
By itself this has little statistical meaning, since our error bars are large. However, our 
distance to SagA* is big enough to yield a sizable effect in the prediction as well as a 
good fit for the average HCC radius. The case of equal radii brings R — 1.7 pc; for 
comparison, the average half-light radius of 157 GCs in the Milky Way is 4.3 pc [40]; 
being mostly heavier, GCs probably arise from agglomerates of HCCs. This value of R 
sets the distance scale of Figs. 1 and 2, and relates the investigated angles to distances 
between 10 and 40 kpc. The HCCs bigger than the diffraction limit appear to be located 
inside the Galactic Halo, in front of the SMC. The less informative LMC map provides 
results roughly consistent with those for the SMC. 

Our SMC-LMC analysis is compatible with a nearly isothermal distribution of 
Herschel cold clouds. If they are distributed in a spherically symmetric fashion 
throughout the Halo, they can explain the full dynamical mass of the Galaxy, provided 
their average mass is of the order of 10,000 M Q . One is thus led to identify the HCCs 
of heretofore unknown origin with the Jeans clusters (JCs) predicted by the Jeans 
instability. Our mass value is in accord with the Jeans instability, which is known 
to have produced globular star clusters, probably from mergers of JCs. The value of 
10,000 Mq is even somewhat smaller than might be expected, which may be due to 
an over-concentration of Jeans clusters in front of the Magellanic clouds. Indeed, some 
clustering, that we observe at small scale, should also be expected on larger scale. 

Thousands of similar cold clouds have been detected by Planck [7] and by 
Herschel/Spitzer [10]. Both the parsec size and 10,000 M weight of the supposed Jeans 
clusters are in accord with the size and gas masses of some of these clouds [8, 9, 10]. 

The amplitude A fixes the ratio Mj C /Rj C , and since in Gravitational Hydrodynamics 
the JCs themselves are isothermal spheres of micro brown dwarfs (/iBDs) of Earth weight 
scale, it determines their velocity dispersion as <x M bd ~ 4 km s _1 , a value known from 
stars in globular star clusters, in support of their presumed connection to Jeans clusters. 

Another surprise is that the deduced cr^bd is close to the speed of sound at 
recombination 5.7 km s _1 , the early moment in the history of the Universe (at 377.000 
years after the Big Bang) when the plasma converted into gas. Intriguingly, the /iBDs 
still seem to move at nearly the speed of sound of the gas that fragmented to create 
them! 

Our analysis points at occultation rather than lensing, which may explain why 
//BDs (MACHOs) can exist but were not observed in microlensing towards the LMC, 
SMC and the Milky Way [41, 42]: they may not act as lenses but as refractive dark 
objects. An occultation event in Centaurus A has lasted 170 days [43]. The cloud's 
mass of 3 - 10 M e agrees with (12); its size of Pluto-orbit scale can be caused by 
heating of the /iBD [16]. A documented doubly-peaked quasar microlensing event [44] 
can be explained as lensing with partial occultation, due to a cloud of about an Earth 
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mass and a size of a few solar radii, located in the Halo. This support our present 
findings, including occultation towards the MCs [45] . 

We conclude that, under the assumption of a spherically symmetric Halo, enough 
cold clouds have been observed towards the Magellanic clouds for the Jeans mechanism 
to explain the full dynamical mass of the Galaxy. In other words, cold JCs can explain 
the dark matter and most of the missing baryons of the Galaxy. The obvious implication 
is that CDM would not be needed to explain dynamical behavior in the Galaxy and, by 
default, at galactic scales. This conclusion was reached before on a the basis of another 
field of the sky and by other arguments [12]. If CDM does not occur in the Galaxy, then 
it likely does not exist at all. But this raises the question: What is the true non-baryonic 
dark matter? A fit to lensing data of the galaxy cluster Abell 1689 has shown that 1.5 
eV neutrinos answer the question [46]. This mass exceeds estimates based on the cold 
dark matter paradigm, but lies below the upper bound of 2 eV. It will be tested in the 
upcoming KATRIN tritium-decay experiment [47] . 

As an outlook, the full data for the SMC will also contain the intensities of the 
HCCs. Investigation of other fields on the sky seems feasible. A dedicated MACHO 
search towards the LMC can look for refractive lensing-occultation events. 
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